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Abstract 



Simplicial quantum gravity has been proposed as a regularization for four dimen- 
sional quantum gravity. The partition function is constructed by performing a 
weighted sum over all triangulations of the 4-sphere. The model is well-defined 
only if the number of such triangulations consisting of N simplexes is exponentially 
bounded. Numerical simulations seem so far to favor such a bound. 
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1. Introduction 



" Simplicial quantum gravity" has been proposed gularization for (Euclidean) 

quantum gravity in four dimensions £[],§■ It was suggested as a natural general- 
ization of similar models studied in two dimensions j|, |3|, [5] and three dimensions 
II 0) H) IE US HH' H11- I n an t nese models the partition functions are given by 



T Or 

where the summation is over all combinatorially inequivalent triangulations of a 
manifold of fixed topology. The factor 1 / CV is the symmetry factor of the triangu- 
lation. It is not important for the following discussion. For pure gravity S(T) is 
bounded by 

\S(T)\<c\T\ (2) 

where |T| denotes the number of simplexes in the triangulation |T|. The inequality 
(0) implies that the partition function is well defined only if the number Af(\T\) of 
triangulations of a fixed topology (which we in the following will assume spherical) 
and consisting of \T\ simplexes is exponentially bounded: 

JV(|T|) < e c|T| . (3) 

In the case of 2d simplicial quantum gravity the exponential bound was proven long 



ago fl4 |. In 3d a number of plausible arguments was given in favor of an exponential 
bound 0, but no complete proof exists so far. Most of the arguments of extend 
to Ad, but again a proof is missing. In the absence of a proof in 3d the question was 
addressed numerically in [0], and good evidence in favor of an exponential bound was 
found. The situation in 4c? seems very similar to the situation in 3d and the finite 
size effects associated with the identification of the smallest possible c in (|3|) were 
even considerably smaller than than in 3d. Recently it has been claimed, however, 
based on new numerical simulations, that the situation in 4d differs profoundly from 
the situation in 2d and 3d and that J\f(\T\) grows factorially |15 |: 



cN(\T\) ~ |T!|V |T| . (4) 
The numerical results they report are in total agreement with earlier simulations 



in the case where we can check (ko = in the notation of [15|, see later) and the 
purpose of this short note is to make clear that there are other, in our opinion more 
plausible, interpretations of the data, which are compatible with the exponential 
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bound (H). In addition we have used the opportunity to extend the simulations to 
somewhat larger volume. 



2. Results 



In 4c? the gravity model ([!]) can be written as |T| 



1 — c -k 4 N 4 (T)+k 2 N2(T) ^ 

T(S 4 ) 



where N±(T) and N 2 (T) denote the number of 4-simplexes and triangles in the 
triangulation T. According to Regge calculus k 2 ~ l/Gc the bare gravitational 
constant. If we have a regular triangulation of S A the number of vertices Nq(T) = 
2 - N 4 (T) + N 2 (T)/2. This allows us to write 



Z ~ £ ^e-^^W no = 2h, n, = h-2k 2 , (6) 



in order to make contact with the notation used in \TE\. If we assume (|3|) there will 
for each value of Kq be a critical value, k1(kq) of k 4 , such that the partition function 
is divergent for k 4 < k%(kq) and convergent for k 4 > k%{kq). The infinite volume 
limit of the system is obtained by letting k 4 approach k 4 (k ) from above. Whether 
this infinite volume limit has an interesting continuum interpretation should then be 
investigated for various values of Kq. A second order phase transition in geometry 
might indicate long range fluctuations and thereby continuum physics. Details of 
how to determine k%{kq) can be found in [0, the only point we want to emphasize 
here is that one determines a pseudo-critical value k2( k o,^4) by analyzing fluctua- 
tions in size around some chosen iV 4 . By changing iV 4 this pseudo-critical point will 
change and in the limit A/4 — > 00 it should converge to k 4 (k ), again provided (0) is 
valid. In case (|J) should be valid there will be no convergence for AT 4 — * 00. 

In order to discuss possible options, let us for simplicity choose k = 0. We are 
then directly testing A/" (A/4). Let us assume that the number of configurations grows 
exponentially. There will be subleading corrections. A reasonable trial ansatz for 
large A/4 is 

log(AT(iV 4 ))~<iV4-giV 4 3 , (3<l. (7) 

This is known to be true in 2d gravity where (3 = (and the power is replaced by 
a logarithm) and (to the extent one can trust numerical simulations) in 3d where 
(3 ~ 2/3. This subleading behavior will reflect itself in the determination of k 4 (A 4 ) 
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which will be given by 

«c(^ 4 ) = K c__£._... j a = i-(3. ( 8 ) 

If we instead assume that (f§) is valid we will get 

Kl(N A ) = a + 6\ogN4. (9) 

Let us first remark that in both CclSCS Ha (A" 4 ) is an increasing function of N4, and in 
case a and 5 are small it might be difficult to distinguish (]|) and @ when only a 
limiting range of A4 is considered. Secondly one should be aware that both (|§|) and 

are asymptotic expressions which might only be valid for large A" 4 . Thirdly, (§) 
and (^D are valid also for kq 7^ 0, the only difference being that k^, a and a might 
depend on Kq, while 5 is independent of kq since ($) implies that various actions 
S(T) can only lead to different exponential corrections. (A4!) 5 is entirely an entropy 
factor coming from the number of triangulations. 

In fig. 1 we show k%{N^ k,q) as presented in [15|] for three values of kq. The curves 



are just linear interpolations between the data points. We see a clear deviation from 
the functional form (|D or (H) for A^ < 4000 in all three cases. Since this deviation 
is qualitatively of the same form in the three cases it is natural to assume that there 
are finite size effects not compatible with (|]) or (|9|) for iV < 4000. In any case it 
seems not reasonable to determine 5 from N4 < 4000 since the functional form here 
shows a strong dependence on k (as also remarked in and 5 is independent 
of kq. Let us now concentrate on the case k = where we can directly compare 
also with our own old (and new) results. In fig. 2 we show K%{Ni) for kq = and 
N± > 4000. When data points can be compared directly there is agreement within 



error bars. The (dashed) straight line is the fit @ found in The fully drawn 
line is a fit to (||) with a = 1/4 (the reason for this choice will be discussed below). 
In both cases there are two free parameters in the fit, a, 5 and k^, c, respectively. 
For the non-trivial parameters, i.e. c and 5 we have c ~ 1.23 and 5 ~ 0.03. 



3. Conclusions 

As is clear from fig. 2 the numerical simulations so far are perfectly consistent with 
an exponential bound on the entropy of triangulations of S 4 . While it is possible 
that one can still make @ compatible with fig. 2 even if we include the new data 
point corresponding to N 4 = 64000, we feel that a fit of the form (||) is more natural 
for > 4000 for several reasons. First it is hard to understand why there should be 
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a factorial growth (N\) s with an extremely small exponent 5 ~ 0.03. As mentioned 
above the origin of such a term is purely combinatorial and if there was a factorial 
growth one would a priori expect 8 ~ 1. Secondly it is difficult for us to understand 
why the situation for S 14 should differ from that of S 3 . Rather, since there is also 
an exponential bound in 2d and 3d it is natural to conjecture that there is an 
exponential bound on the number of triangulations of S d for any d. 

It is interesting to note that data are compatible^ with the choice a = 1/4. This 
means that the exponent (3 in (|7|) is 3/4 and if we tentatively write N4 = L 4 eq. (0) 
reads: 

logA/"(L 4 ) = K4L 4 — cL 3 . (10) 

In 3d we found |7j (3 ~ 2/3 and introducing in 3d the notation N 3 = L 3 one has the 
analogue of (|T0|): 

logA/"(L 3 ) = - cL 2 . (11) 

This kind of behavior might be related to the fact that the generic triangulation 
for /?o = 0, both in 3d and 4d is very singular. There are vertices of very high 
order (almost comparable with |T|) and the d — 1 dimensional boundary of the 
unit ball surrounding such a vertex could maybe act as "pseudo-boundary" of the 



closed (i-dimensional manifold. It would explain the boundary-like terms in (|l0l ) 
and (|TTD. When k,q is increased (decreasing bare gravitational constant) the vertices 
of very high order disappear. At the same time the value of /3 seems to decrease, 
maybe approaching zero at the observed phase transition^. In addition it would also 
explain why the corresponding term cL is absent in 2d. The distribution of high- 
order vertices in 2d can be calculated analytically and the probability of having such 
vertices falls off exponentially with the order of the vertex. 



2 We should note that a is not very well determined from a fit to the data shown in fig. 2. Any 
a in the range 0.2 < a < 0.3 will give a reasonable fit. 

3 The value of j3 (or a) is not very well determined from the fits to k%{N4) by the present data 
for kq near or above the phase transition point. In 3d the situation is better and it seems as if 
f3 w at the transition |l3| . The study of baby-universe distributions also supports this 
picture. 
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Figure caption 



Fig. 1 The data points taken from ||15||. The plot is log N 4 versus k^N^, k ) for k 



(lower curve), k = 0.5 (middle curve) and k = 1.0 (upper curve). 

Fig. 2 log A^ 4 versus k%{N±, k ) for k = 0. Circles indicate the data points from the 
plot in fl5| ) the size of the circles being approximately the error bars. The 
dots (with error bars) indicate our own data. We see that there is perfect 



agreement. The dashed straight line is the fit to (g) as presented in ||15|| . The 
fully drawn line is the fit to (||) with a = 1/4. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-lat/9405010vl 



